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Abstract 


The dwarf planet Pluto is known to host an extended system of hve co-planar 
satellites. Previous studies have explored the formation and evolution of the sys¬ 
tem in isolation, neglecting perturbative effects by the Sun. Here we show that 
secular evolution due to the Sun can strongly affect the evolution of outer satel¬ 
lites and rings in the system, if such exist. Although precession due to extended 
gravitational potential from the inner Pluto-Charon binary quench such secular 
evolution up to acrit ~ 0.0035 AU (~ 0.09 Rmii the Hill radius; including all of 
the currently known satellites), outer orbits can be signihcantly altered. In par¬ 
ticular, we hnd that co-planar rings and satellites should not exist beyond acrit] 
rather, satellites and dust particles in these regions secularly evolve on timescales 
ranging between 10^ — 10® yrs, and quasi-periodically change their inclinations 
and eccentricities through secular evolution (Lidov-Kozai oscillations). Such os¬ 
cillations can lead to high inclinations and eccentricities, constraining the range 
where such satellites (and dust particles) can exist without crossing the orbits 
of the inner satellites, or crossing the outer Hill stability range. Outer satellites, 
if such exist are therefore likely to be irregular satellites, with orbits limited to 
be non-circular and/or highly inclined. These could be potentially detected and 
probed by the New-Horizon mission, possibly providing direct evidence for the 
secular evolution of the Pluto satellite system, and shedding new light on its 
origins. 
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Introduction 


Since its discovery in 1930 flTombaughl Il946l) the dwarf planet Pluto has provided us 
with valuable information on the origin of the Solar system, its structure, its dynamical 
evolution and its basic building blocks. It was the hrst object to be discovered in the Kuiper 


belt and suggest its existence; its orbital resonan ce with b 


provided strong evidence for Neptune migration (iMalhotra 


eptun e and high inclination 


I 993 I) , and thereby provided 


important clues into the early hi story of the Solar System an d its assemblage. The discovery 


of its close companion Charon fjChristv &: Harrington 


19781) revealed the existence of a 


new type of planetary systems, binary planetesimals/dwarf-planets, and i 


tidally-locked conhguration provided further clues for its collisional origin fjCanup 


s close and 


209511 . 


A series of discoveries in recent years showed the existence of no-less than four 
additional satellites, including Styx, Nix, Kerberos and Hydra orbiting Pluto on circular 
orbits with semi-major axes ranging between ~42K-65K km, i.e. residing between 
~ 0.007 — 0.011 of the system Hill radius, Rmii = ap(l — ep){mp / SrriQy^^ ; where ap, ep 
and mp are the semi-major axis (SMA) of the orbit of Pluto around the Su n; the orbital 


eccen tricity; and the mass of Pluto, respectively; tuq is the mass of the Sun flBrozovic et ah 


2015 


the 


2006 


o rmation, stability and evolution o 


Stern et ah 


2006 


Ward fc CanuD 


this e xtended system 


2006 


Lithwick fc Wu 


to a series of studies on 

{Brown 

2002 

Lee & ’ 

^eale 

- 

00 

0 

0 

CM 

b 

CanuD 

2011 



Youdin et al.l 2012: Giuliatti Winter et ah 2013: Cheng et al.ll2014bl: Giuliatti Winter et al 


2014: 


Kenvon fc Bromlev 


Porter fc Stern 


2014 : 


Bromlev &: Kenvon 


2015 


Giuliatti Winter et al. 


2015 


20151 ). that could shed light not only on the origin of the Pluto system. 


but also provide new clues to the understanding of the growth of planetary system in 
general, as well as on migration and resonant capture processes in planetesimal disks. 
Better understanding of the properties of the Pluto system and its constituents is therefore 



















































































4 


invaluable for advancing our knowledge of the Solar System, its building blocks and and 
its origins. The New-Horizon mission, launched in 20 06, was designed t o explore the 


Pluto-Charon system and help accomplish these goals flYoung et al 


20081). New horizon is 


now reaching its climax, as it makes its close flyby near the Pluto-Charon system, aiming 
to characterize Pluto and its satellites at an unprecedented level, and potentially detect 
additional satellites and/or planetesimal rings in this system. Understanding the dynamical 
history and the current configuration of the Pluto satellite system system are therefore 
invaluable for realizing New-Horizons’ data collecting and characterization potential, and 
use them in the interpretation of the data and their implications for the origins of Pluto 
and the Solar System history. 


Here we show that secular evolutionary processes due to gravitational perturbations by 
the Sun, previously neglected, may have played a major role in the initial formation of the 
Pluto system as well as in the later evolution of its satellite system and its conhguration 
to this day. As we discuss below, our hnding provide direct predictions for the possible 
orbits of outer moons in the Pluto-Charon system, if such exist. In particular we show 
that co-planar planetesimal rings can not exist beyond a specific critical separation from 
Pluto-Charon (significantly smaller than the Hill radius of the system), and that any moon 
residing beyond this critical separation must be an irregular, inclined and/or eccentric moon. 
Moreover, we constrain the orbital phase space regimes (semi-major axis, eccentricity, 
inclination) in which dust, planetesimal-rings and/or outer moons can exist on dynamically 
and secularly stable orbits. 


2. Secular evolution of sattelites in the Pluto-Charon system 

Charon is Pluto’s closest and largest moon, with semi-major axis of ape ~ 19.4iP km 
(period of 6.387 days) and mass ~ 1.52 x lO^^g (mass ratio of fi = rric/mp = 0.116). It 
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revolves around Pluto on a tidally-locked (i.e. in double synchronous state) circular orbit 
with an inclination of / ~ 119° in respect to the orbit of Pluto around the Sun (which 
in itself is inclined in an angle of ~ 17° in respect to the ecliptic plane). Together with 
the Sun, the Pluto-Charon system comprises a hierarchical triple system, with Pluto and 
Charon orbiting each other in an inner binary orbit, and their center of mass orbiting the 
Sun in an outer binary orbit. The evolution of such hierarchical gravitating triple systems 
may give rise to secular dynamical processes, in which the outer third body perturbs 
the orbit of the inner binary, leading to secular precession of the inner binary, typically 
exciting its eccentrcity. In particular, triples with high mutual inclinations between the 
inner and outer orbits (typically >~ 40°, and somewhat lower for eccentric syste ms) 
are subject to qua si-periodic secular evolution, so called Lidov-Kozai (LK) cycles flLidov 


1962 


Kozai 


196211 . w 


binar y planetesimals 


rich could affect Plu to’s satellites as well as other S olar System 


Perets &: Naoz 


20091) and satellites of the gas-giants fjCarruba et ah 


20021). Such LK evolution leads to the inner orbit precession and in turn gives rise to 


high amplitude eccentricity and inclination oscillations that occur on the LK precession 
timescales, which are much longer than the orbital period. 


In principle, the high mutual inclination of the Charon-Pluto system makes it highly 
susceptible to such LK secular evolution, and in the absence of other forces it would have 
lead to a collision between Pluto and Charon. However, such processes are highy sensitive to 
any additional perturbations of the system precession. In particular, the effects of the tidal 
forces between Pluto and Charon at their current separation lead to a significant precession 
of their orbital periapse, which will be denoted gpc- This, in turn, quenches any LK secular 
evolution, and keeps Pluto and Charon on a stable circular orbit, with no variations of the 
orbital eccentricity and the mutual inclination in respect to the orbit around the Sun. 


The strong dependence of the tidal forces on the satellite separation renders the tidal 
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effects on Pluto’s four other known satellites negligible. However, periapse precession can 
still be induced by the non - point-mass gravita tional-potential of the inner Pluto-Charon 


binary ( 


^ee &: Peak 


2006 


planets (iHaniers et ah 


H 


2015bl: 


amers et al. 


2015al). similar to the case of circumbinary 


Martin et al 


12015 


Munoz&Lai 


2015fl . In this case, if the 


inner-binary induced precession timescale, rpc, is shorter than the LK precession timescale, 
tlk, the orbit will not librate and the LK secular evolution will be quenched. If, on 
the other hand, LK period is shorter than the precession timescale, the binary-induced 
precession is slow and the orbit has sufficient time to librate and build-up signihcant LK 
oscillations. The critical SMA at which these timescales become comparable is given by 
(Fig. 1; see appendix for full derivation) 


^crit — 


3/2 , , n 1/5 


3 ap^apTUpmc (l - e^) (50^ - 1) 

8 (mp-f-me) (1 - e^)^mo 


( 1 ) 


where parameters with sub-indexes PC correspond to the orbital parameters of the 
Pluto-Charon orbit, and the m index refers to the orbital parameters of the orbiting moon. 
The inner orbit semi-major axis (SMA) for Pluto-Charon mutual orbit is opc; the SMA of 
the moon orbit around Pluto-Charon is given by am] the eccentricity of the inner orbit, is 
epc] the eccentricity of the outer orbit is e^; the arguments of the pericenter of the inner 
and outer orbits are gpc, Qm, respectively and 9 = cosi is the cosine of the inclination 
between the two orbit planes denoted by i . Charon mass is given by me and the mass of 
the Sun is denoted by m©. 


In Fig. [T] we compare the satellite precession timescales due to the inner-binary 
induced precession with the precession timescales due to LK secular evolution. The opposite 
dependence of these two precession timescales on the orbital separation give rise to a critical 
separation at which they equalize, as derived in Eq. [H Beyond the critical separation,Ocrit, 
at which the timescales become comparable, satellites/rings become highly susceptible to 
LK secular evolution. The effect of the inner binary can be combined with the full secular 
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Fig. 1.— Comparison between the precession timescales for a satellite orbit induced by perturbations by the 
Sun (Lidov-Kozai evolution) and the precession induced by the inner Pluto-Chaon binary. Beyond the critical 


separation,Ocrit, at which the timescales become comparable, satellites/rings become highly susceptible to 


Lidov-Kozai secular evolution. 












Fig. 2.— Examples for the secular evolution of outer satellites in the excluded regions. Top: The evolution 
of a satellite on an initially co-planar (with respect to the Pluto-Charon orbit;119° with respect to the orbit of 
Pluto around the Sun), circular orbit. The eccentricity and inclination of the orbits change periodically due 
to LK evolution until the satellite separation becomes too large and crosses the Hill stability region, making 
the satellite orbit unstable. Right: The evolution of a satellite on an initially inclined (-20 degrees in respect 
to Pluto-Charon; 99° with respect to the orbit of Pluto around the Sun), circular orbit. The eccentricity and 
inclination of the orbits change periodically due to LK evolution until the satellite separation becomes too 
small and crosses the orbit of Hydra; at which point the satellite system is likely to destabilize and/or the 
satellite can collide with Hydra; such satellite are therefore unlikely to exist in the Pluto-Charon system. 
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Fig. 3.— Allowed and excluded regions for the existence of outer satellites and dust particles. The minimal 
and maximal separations of satellites during their secular evolution are shown as a function of their initial 
separation from the center of mass of the Pluto-Charon system. Top left: Satellites on initially co-planar 
circular orbits do not not secularly evolve up to distances comparable with the critical separation. Beyond 
this point the eccentricities and inclinations of such satellites are periodically highly excited. The secular 
evolution drives the satellites through a range of separations, enclosed by the blue solid lines. Satellites at 
sufficiently high eccentricities become unstable as they cross the stability region (Hill stability; ^ 0.5Rh for 
retorgrade orbits; top dashed vertical line) at apocenter, or if they cross the orbits of the inner satellites at 
pericenter (bottom dashed line and left vertical lines show Hydra SMA/separation). Regions in which such 
instabilities occur are therefore excluded regions where outer satellites can not survive. Top right: Similar 
but for satellites on initially inclined orbits of +10°. Bottom left: The same for satellites on initially inclined 
(-20 degrees with respect to Pluto-Charon orbit; 99° with respect to the orbit of Pluto around the Sun) 
circular orbits. Bottom right: The same for retrograde orbits (inital inclination of —55°). 
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equations of motion of the quadruple system (Pluto-Charon+outer satellite+the Sun) to 
derive the full orbital evolution of the outer satellite in time. We use such derivations to 
hnd the maximal inclination and eccentricities attained by an initially co-planar satellite 
orbiting Pluto-Charon on a circular orbit, as a function of the satellite initial separation 
from Pluto (Fig. [2]). As can be seen, the LK precession timescale is longer than the 
binary-induced precession timescale for satellites at small separations, including the current 
locations of Pluto satellites. LK-evolution induced eccentricity and inclination oscillations 
are completely quenched in these regions below acrit, allowing for the stable orbits of the 
currently known satellites. However, beyond this point, LK oscillations are not (or only 
partially) quenched. Satellites in these regions experience large amplitude eccentricity and 
inclination oscillations. In principle, satellites excited to high eccentricities might even cross 
the orbits of the inner satellites at peri-center approach or extend beyond the Hill stable 
region at apocenter. The former case leads to strong interactions with the inner satellites; 
such satellite will eventually destabilize the satellite system or collide with one of the inner 
moons. In the latter case the perturbations by the Sun will destabilize the orbit of the 
moon, likely ejecting it from the system or sending into crossing orbits with the inner moons 
(see examples in Fig. [2]). Regions in which such instabilities occur are therefore excluded 
regions where outer satellites can not survive (Fig. |3]). 

We can therefore exclude the existence of moons with some given orbital parameters 
a, i, e li they do not follow the stability criteria 

(^hydra ^ ®[l L c)] andci[l -(- Gmaxifl) L c)] < Rhill 

where emax{(^,h^) is the maximal eccentricity of the moon during an LK cycle (which can 
be analytically or numerically derived from the coupled secular evolution equations of the 
orbital parameters; including both binary precession and LK precession terms; see Fig. 
H^nd appendix). Note that not far beyond the critical separation the LK timescale becomes 
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much shorter than the binary-precession timescale, at which point the secnlar evolntion is 
completely dominated by LK evolntion, and the binary precession terms can be neglected. 


i.e. for snch large separations one can nse t 


le an alytic solntion for Cmax dnring LK 


oscillations, Cmax = [1 “ (5/3) cos^/]^/^ (ILidov 


19621) which is independent of the separation. 


This criterion therefore provides a robnst map of the exclnded and allowed orbital phase 
regions in which onter moons may exist in the Plnto-Charon system. Moreover, the cnrrent 
stable conhguration of the known inner satellite system can therefore be nsed to constrain 
the region where farther satellites can exist. 


For initially circular co-planar moons with a > acru, ^max — 0.78 and we get 
acrit(l — ^max) > cihydra, i-G. such moous uever evolve to hydra crossing orbits. This can 
be seen in Fig. [3] showing that initially regular moons on co-planar circular orbits beyond 
am_crit evolve through LK cycles and can obtain high inclinations and eccentricities, but 
their trajectories never cross the orbits of the inner moons. Such moons, even if they 
were initially formed on co-planar circular orbits (“regular” moons) can not sustain such 
orbits; they secularly evolve, and if observed they are expected to be irregular moons, 
with non-negligible eccentricity and/or inclination in respect to the Plnto-Charon orbit. 
Moreover, initially regular outer moons can exist only in a limited outer region beyond the 
critical separation, as satellites at even larger separations cross the Hill stability radius. 
Large regions of the orbital phase space are therefore excluded by these criteria for the 
existence of moons. 


Our results therefore constrain the maximal extension where regular moons and 
planetesimal disks can exist around the Plnto-Charon system. We can therefore predict that 
any moons residing beyond the critical separation, if such exist, will be eccentric/inclined 
irregular moons, and we map the specihc orbits allowed for such moons. These hnding 
also show that secular LK evolution can not be neglected in studies of the formation and 
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semi-major axis [AU] 


Fig. 4.— The maximal eccentricity attained by an outer moon as a function of its semi¬ 
major axis (for satellites on initially co-planar circular orbits). The maximal eccentricity 
excited during the secular evolution depends on the the coupled effect of the perturbations 
by the Sun (LK evolution) and the precession induced by the Pluto-Charon inner binary. In 
the inner regions LK evolution is quenched and the satellite keeps its initial eccentricity and 
inclination, at the outer regions binary precession become negligible and the maximal eccen¬ 
tricity is derived directly from the LK evolution. In the intermediate regimes both process 
are important and the maximal eccentricity can be derived numerically or be approximated 
analytically. 









evolution of the Pluto-Charon system and its satellites, as hitherto (tacitly) assumed. 



- 14 - 

REFERENCES 

Blaes, O., Lee, M. H., & Socrates, A. 2002, ApJ, 578, 775 

Bromley, B. C. & Kenyon, S. J. 2015, arXiv: 1503.06805 

Brown, M. E. 2002, Annnal Review of Earth and Planetary Sciences, 30, 307 

Brozovic, M., Showalter, M. R., Jacobson, R. A., & Bnie, M. W. 2015, Icarus, 246, 317 

Canup, R. M. 2005, Science, 307, 546 

—. 2011, AJ, 141, 35 

Carruba, V., Burns, J. A., Nicholson, P. D., & Gladman, B. J. 2002, Icarus, 158, 434 

Cheng, W. H., Lee, M. H., & Peale, S. J. 2014a, Icarus, 233, 242 

Cheng, W. H., Peale, S. J., & Lee, M. H. 2014b, Icarus, 241, 180 

Christy, J. W. & Harrington, R. S. 1978, AJ, 83, 1005 

Ford, E. B., Kozinsky, B., & Rasio, F. A. 2000, ApJ, 535, 385 

Ciuliatti Winter, S. M., Winter, O. C., Vieira Neto, E., & Sfair, R. 2013, MNRAS, 430, 
1892 

—. 2014, MNRAS, 439, 3300 
—. 2015, Icarus, 246, 339 

Hamers, A. S., Perets, H. B., Antonini, F., & Portegies Zwart, S. F. 2015a, MNRAS, 449, 
4221 

Hamers, A. S., Perets, H. B., & Portegies Zwart, S. F. 2015b, arXiv:1506.02039 


15 


Kenyon, S. J. & Bromley, B. C. 2014, AJ, 147, 8 

Kozai, Y. 1962, AJ, 67, 591 

Lee, M. H. & Peale, S. J. 2006, Icarus, 184, 573 

Lidov, M. L. 1962, Planet. Space Sci., 9, 719 

Lithwick, Y. & Wu, Y. 2008a, ArXiv0802.2951 

—. 2008b, arXiv:0802.2939 

Malhotra, R. 1993, Nature, 365, 819 

Martin, D. V., Mazeh, T., & Fabrycky, D. C. 2015, ArXivl505.05749 
Michaely, E. & Perets, H. B. 2014, ApJ, 794, 122 
Munoz, D. J. & Lai, D. 2015, arXiv:1505.05514 

Naoz, S., Farr, W. M., Lithwick, Y., Rasio, F. A., & Teyssandier, J. 2013, MNRAS, 431, 
2155 

Perets, H. B. & Naoz, S. 2009, ApJ, 699, L17 
Porter, S. B. & Stern, S. A. 2015, arXiv: 1505.05933 

Stern, S. A., Weaver, H. A., Steffi, A. J., Mutchler, M. J., Merline, W. J., Buie, M. W., 
Young, E. F., Young, L. A., & Spencer, J. R. 2006, Nature, 439, 946 

Tombaugh, C. W. 1946, Leaflet of the Astronomical Society of the Pacific, 5, 73 

Valtonen, M. & Karttunen, H. 2006, The Three-Body Problem 


Ward, W. R. & Canup, R. M. 2006, Science, 313, 1107 


- 16 - 

Youdin, A. N., Kratter, K. M., & Kenyon, S. J. 2012, ApJ, 755, 17 

Yonng, L. A., Stern, S. A., Weaver, H. A., Bagenal, F., Binzel, R. P., Buratti, B., Cheng, 
A. F., Cruikshank, D., Gladstone, G. R., Grnndy, W. M., Hinson, D. P., Horanyi, 
M., Jennings, D. E., Linscott, I. R., McComas, D. J., McKinnon, W. B., McNntt, 
R., Moore, J. M., Mnrchie, S., Olkin, C. B., Porco, C. C., Reitsema, H., Renter, 
D. C., Spencer, J. R., Slater, D. C., Strobel, D., Snmmers, M. E., & Tyler, G. L. 
2008, Space Sci. Rev., 140, 93 


Appendix 


A. Calculation of acrit 

In the canonical two-body problem the bodies are treated as point particles. Any 
deviation from the point mass approximation changes the solution, and the orbital 


evolution. For example, if we change the primary from a point mass to an o 


the resulting orbit of the secondary will be a precessing ellipse flCheng et al 


jlate s phere 


2014ah . The 


same end result will be the case if we were to change the primary point mass to a stable 
binary system. In this case the secondary orbits the center of mass of the inner binary 


in a precessing Kep 


erian ellipse, as is the case for Pluto-Charon and their coplanar zero 


eccentricity moons (ILee fc Peald 120061) . In the followi ng we describe the precession rate 


the more genera l case. Following similar c 


Naoz et al. 


2013 


Michaelv &: Perets 


20 lJ; 


erivations flFord et al. 


Hamers et al 


20001 


3 precession 

rate 

Blaes et al. 

2002 


or 


2015al) we find that the precession 


rate of the tertiary in the quadrupole expansion, the outer orbit, due to the inner orbit, is 


This manuscript was prepared with the AAS IATeK macros v5.2. 
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given by (the inner Pluto-Charon orbit is constant): 


= 3(^2 1[2 + elc (3-5 cos 2gpc)\ 


dt 


G 


PC 


where 


+ 3^2 <* — [4 + Qelc + (50' - 3) (2 + 3elc - cos 2gpc)] 


Co = 


Gmpmcmr, 


16 (mp + me) am (1 - 
the inner and outer angular momenta are 


apc 


(Al) 


(A2) 


Gpc — 


mpmc 
mp + me 


\JG {mp + me) ape (1 - e|,^) 


(A3) 


mm {mp + me) 7777-^^^- T. -yr 

Gm = -^^- \/G [mp + me + mm) am (1 - e^) 

mp + me + m^ 

For Pluto-Charon case epc = 0, due to tidal interaction, therefore eq. reduces to 

d^r 


(A4) 


dt 


= 6^2 


G 


29 502 - 1 

+ 


PC 


Gr 


(A5) 


Equation (IA5I) sets a timescale for the orbit precession. This timescale quenches Kozai 
cycle acting upon the moon. In this case Kozai mechanism acting upon the following 
triple system, the inner binary is the moon orbiting Pluto-Charon. In this approximation 
Pluto-Charon considered to be a single object, due to tidal interaction that sets their 
orbital parameters. The outer orbit is the orbit that Pluto orbit the Sun. In order to hnd 
the SMA on which the timescale are equal and hence hnd the critical SMA from which the 
moon orbits is governed by Kozai mechanism we need to equate the Kozai timescale and 
the precession timescale and solve for am, the inner orbit SMA: 

f 20 , 502 - 1\ fdgr, 


6Go 




Gpe 


+ 


Gr, 


dt 


(A6) 


Kozai 
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where {dgm/dt)Kozai is given by 

\ dt ) (jYip rnc + (1 - 

where Cp is Pluto’s orbital eccentricity and G is Newton’s constant. In the general case we 
can solve this equation numerically. However, in the test particle limit we can neglect the 
first term in the parenthesis of the left hand side and solve for am analytically. We can see 
this fact if we compare the importance of the two term: 


2e/Gpc _ f ‘^0 \ 

{5e^-i)/Gm 


(A8) 


where k is constant with units. We can see the in the test particle limit '—> 0 the hrst 

term in negligible except an extreme case. In the case of a specific moon inclination that 
satisfies the following condition 

502 - 1 = 0 (A9) 


which corresponds to 


cos^ i 


1 

5‘ 


(AlO) 


If the approximation is fulfilled then we can solve for and get the critical 

3 a?p(ja\,mpmc (l - (502 - 1) 

8 (mp+ mc) (1 - e^)^mo 


(All) 


Close inspection of (lAllD we notice an interesting feature of the dynamics. Due to the term 
(1 — e^) in the denominator moon orbits with high orbit eccentricity are more stable, with 


respect to Kozai oscillations than circular ones, which is counterintuitive. 


B. Calculation of maximal eccentricity for initially co-planar circnlar orbits 

The maximal eccentricity obtained by a moon during its long term evolution can be 
derived by solving the full coupled secular evolution equations. However, one can derive 
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approximate analytic solutions in some specific cases. In the following we derive the 
maximal eccentricity for a moon on an initially co-planar circular orbit around the inner 
binary. 


In order to find the maximal eccentricity for initial nearly circular orbit of the moon’s 
orbit around Pluto-Charon we write the precession equation of motion of the moon’s orbit 
explicitly. The LK precession rate in the quadrupole expansion in the nearly circular initial 
orbit can be added to the binary precession by the inner Pluto-Charon system to give 


dQm 

dt 


dQm 

dt 


+ 


percession 


d(]m 

dt 


(Bl) 


Kozai 


dQm 

dt 


3 G^/‘^mpmcapQ (50^ -1) 3 G^l'^msam /„ p. . 2 -20 

o- o/n 7/2 -^ + 7 - 77 -(2 - 5sm ^f^sin %) 

^{mp^mcf^ arL (1 - e^) ^{mp + mc) ^ al{l-e],f^ 

(B2) 


and from the definition of acrU) see equation flAlip we can substitute and get 


dQm 

dt 


(S/2)n 

Txozai 


(2-5 sin^ Qm sin^ '^) + t 

' 3 \ am J 


if we redehne t into a more convenient coordinate 


T = t ■ 


(3/2) .r 

Txozai 


(B3) 


(B4) 


(B5) 

an immediate result is that as am becomes larger than acrit the precession rate tends to the 
standard LK rate and therefore the standard Cmax- The standard LK treatment for the test 
particle problem leads to an additional constant of motion 


dQr. 


dr 


= 2 — 5 sim Q, 


. 2 • I ^ \ 

m sm % + - ' - ' 


3 \ am J 


2 = 5 sin^ Qm sin^ i 


(B6) 
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(for derivation see, e.g. IValtonen fc Karttunenll2006l ). In our system the correction to the 
maximal eccentricity is due to the term that proportional to am_crit/am- Therefore the 
constant of motion is 

2 + = Ssin^^fmSin^i (B7) 

o 

where we dehne q = acrit/(^m- From the former equation we can get the relation between 


Qm and sinh 
see Ford et ah 


bgether with the secula r evolution of the orbit eccentricity (for full equation 


2000 


Blaes et ah 


2 OO 2 I) we get the following condition 


. 2 . 2 + 

sm t > - 


(B8) 


and by using the conservation of inner orbit angular momentum equation 


1 - el cos *0 = V^-elnax\ 1 - 




2 + |g 


cos^zo = (1 


3-|g5 


(B9) 

(BIO) 


1 — COS^ Zo 




and therefor the maximal eccentricity is 


Gmax = -I / 1 - COs2 Zq 




(Bll) 


(B12) 


note that for 




(B13) 


the maximal eccentricity is ill defined; in these cases is well below acrit and the secular 
LK evolution is completely quenched. 


























